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Physics potential of the decays B s d — » and B s — » J/if/<f> in the ATLAS 
experiment at the LHC 
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The CP asymmetry predicted by the Standard Model in the decay modes B s — » J /if/ + 77 and B s — > J/i// + <f> is very small, and the obser- 
vation of a sizeable effect would be a clear indication of new physics beyond the Standard Model. Contrary to the final state J/if/Q, the 
former final state is a CP eigenstate, and therefore complex angular distribution analyses are not necessary for the asymmetry measure- 
ment. In this paper, the approaches for the study of these two decay channels are presented. 



1 Introduction 

The decay mode B s — > J/if/Tj is analogous to the mode 
B s — > J/if/(f>, which has been studied extensively in view of 
CP violation measurements. In these two decay modes, the 
CP asymmetry predicted by the Standard Model is very 
small, and the observation of a sizeable effect would be a 
clear signal of physics beyond the Standard Model. 

The decay channel B s — > i/if/rj can be used to measure var- 
ious parameters in the B s -meson system. In the ATLAS 
experiment, assuming an integrated luminosity of 30 fb _1 , 
roughly 10 000 decays can be reconstructed with a signal- 
to-background ratio of about 1:1, as will be shown later in 
this paper. It can be thus foreseen that the measurement 
of the branching fraction of this mode, the measurement of 
the B s lifetime, as well as other measurements will be fea- 
sible, which constitues an important cross-check of other 
measurements performed using other decay channels. 

Moreover, the decay mode B s — > i/if/rj, when combined 
with its U-spin symmetric channel, Bd — » i /<///], provides a 
new strategy for the extraction of the y angle of the Unitar- 
ity Triangle [[Q. Indeed, while the phase exp(/y) is CKM 
suppressed for the decay mode B s — > 3/if/rj, this is not the 
case for the decay mode Bd — > 3/if/rj. CP violation effects 
could thus potentially be more easily accessible through 
the latter mode by measuring the time-dependent asymme- 
try of Bd. The overall normalization can be fixed by mea- 
suring CP averaged rates of Bd — > J/if/J] and B s — » i/if/rj, 
which assumes the validity of the U-spin symmetry (the 
approximate symmetry of u, d and s quarks). The measure- 
ments seem to be, however, out of reach due to the small 
expected branching ratio of Bd — > i/if/rj, and the overlap- 
ping mass peaks of the B s and Bd rl2l. 

The B s - S s system is characterized by two eigenstates 
with different masses and decay rates. The experimental 
determination of the mass difference Am s and the rate dif- 
ference AT S will be valuable input for flavour dynamics in 
both the Standard Model and its possible extensions. While 
the measurement of Am s has been proven to be accessi- 



ble through the decays B s — > D s 7r and B s — > D s ai, AT S is 
expected to be precisely measured in ATLAS through the 
decay channel B s — > J/if/<p. 

2 The B s d — > J/if/rj decay channel 

Signal samples of B S d — » were generated using the 
Monte Carlo program PYTHIA 5.7 [0. The J/i/r meson 
was forced to decay into p + p~, and only events passing 
the ATLAS level- 1 trigger requirements for B hadrons [ 
4 1 (a muon with a pj > 6 GeV and \rj\ < 2.4) 1 and the 
77 meson decaying into yy were retained. The branching 
ratio for the decay rj — » yy is (39.33 + 0.25) %[ 5 1. Events 
were further selected to satisfy the ATLAS second level 
trigger: the presence of a second muon with pj > 3 GeV 
and I77I < 2.5. The resulting signal samples were consisting 
of about 17,000 B s events and 15,000 B d events. 

A full GEANT-based simulation was used to simulate the 
response of the ATLAS inner detector and the electromag- 
netic calorimeter. The muon chambers were not included 
in the simulation, but only real muons were included in 
the analysis. The muon identification efficiencies were 
assumed to be 85% for the first muon (the muon with 
pi > 6 GeV) and 78% for the second muon (the muon 
with pj > 3 GeV). The 85% efficiency for the muon with 
the higher transverse momentum includes both reconstruc- 
tion and level- 1 trigger efficiencies. 

The J/if/ was reconstructed by fitting pairs of opposite- 
charge muons passing the level-2 selection (p6p3) 2 to a 
common vertex and by calculating their invariant mass. 
The^ 2 /d.o.f. of the vertex fit was required to be less than 
6.0, and the decay vertex of the J /if/ was required to be de- 
tached from the primary vertex by at least 250 pm in the 
transverse plane. The J /if/ mass resolution was found to be 

1 Throughout this paper, the symbol pj is used for the transverse momen- 
tum with respect to the beam direction, and )/ for the pseudorapidity. 
2 In the following, the notation fi6 (fi3) is used to indicate a muon with a 
px > 6 (3) GeV. 



2 



Workshop on the CKM Unitarity Triangle, IPPP Durham, April 2003 




T| Invariant Mass (GeV) 



Figure 1. Invariant mass distribution of ?/-mesons in signal events 
after all the cuts. 

39 MeV. The resolution of the J/i// decay vertex in the trans- 
verse plane was 64 pm [|SJ. The efficiency of the recon- 
struction within a three-standard-deviation mass window 
was 79%. If a J/0- was found, the analysis proceeded with 
77-reconstmction in the electromagnetic calorimeter. The 77 
reconstruction is described in more detail in Ref. [ 7 1. The 
mass distribution, Fig. 1, was fitted with two Gaussians, 
and the obtained mass resolution was cr -70 MeV. The 
total reconstruction efficiecny was 2.3% within two stan- 
dard deviations from the nominal mass. 

The }/iff and rj candidates were combined to reconstruct 
both the B s and Bj. A mass resolution of around 67 MeV 
was found for both channels, with an overall signal effi- 
ciency of ~ 1 % within three standard deviations from the 
nominal mass. 

The background considered was a sample of B — > J/tfrX 
decays, processed and selected in the same way as the sig- 
nal events. The resulting efficiency was less than 0.009%. 
The table below summarizes the event rates for both sig- 
nals and background. 

For the measurement of the CP asymmetry in the B s - 
system, the so-called Same-side jet-charge tagging method 
was used for distinguishing production of B s and B s rl6l. 

The jet charge was defined as : 

^i el ~ "v F' 
Li Pi 

where q, is the charge of the 2 th particle in the jet, and p t is 
the momentum. All charged particles with pj > 0.5 GeV, 
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Figure 2. Invariant mass distribution of the reconstructed B s - 
signal (white) and the background (dashed). 





Op = -10° 


e P = -20° 


Sr(B s -» 

jyobs 


8.3 ■ 10~ 4 
8 400 


9.5 • 10~ 4 
9 600 


Sr(B d -» J/^;) 

.yobs 


4.1 • 10" 6 
200 


1.6 - 10" 6 
80 


jyobs 
back 


10 800 


10 800 



Table 1. Assumed branching ratios for B s , t | — » J/if/T], and the esti- 
mated numbers of signal and background events for an integrated 
luminosity of 30 fb . 8 P is the rj — rf mixing angle ( see Ref. [ 

ED)- 

\rj\ < 2.5 around the reconstructed B-meson were included 
in the jet, if the distance between the particle and the B- 
meson, AR = ^jrj 2 + (f> 2 , was less than 0.8. Particles not 
originating from near the primary vertex were excluded 
by requiring that the transverse impact parameter do was 
less than 1 cm and the z-coordinate of the particle at the 
point of closest approach was within 5 cm of the z of the 
primary vertex. The particles from the B-meson decay 
itself were excluded as well. In the analysis, the recon- 
structed B-meson was defined as B s (B s ) if the jet-charge 
had Qjet > + c (Gjet < where c is a tunable cut. 

The exponent k was optimized to maximize the tagging 
quality factor Q = D^ ag ■ e Vdg , where D tag is the tagging 
dilution factor, originating from wrong-sign tags (A ag = 
1 - 2W, where W is the fraction of wrong-sign tags), and 
e tag is the tagging efficiency. The optimum parameters were 
found to be k — 0.2, c - 0.2, resulting in a quality factor of 
3.85% with Aag = 0.26, e tag = 0.57. 
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The observable asymmetry is : 

flobs(0 = Dacp(t) = D sin <pM sin Am s t, 

where D = Aag-Dback combines the experimental dilu- 
tion factors due to mistagging and background. Db ac k = 
N^ bs /(N^ bs + N^ k ) is the ratio of the number of observed 
signal events to the total number of observed events. 

Assuming that the decay time resolution is At =0.073 ps, 
the error on the CP asymmetry was estimated to be: 



£(sin 4>m) 
£(sin 4>m) 



0.27 ', x s — Am s /T s 
0.31, x 5 = 30, 



19, 



where 
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where A = sin#c, 8c is the Cabibbo angle, tj is the height 
of the Unitarity Triangle and y is one of the angles of the 
Unitarity Triangle. 



Figure 3. <j> s - x s regions allowed by the Standard Model predic- 
tions and predictions from Left-Right Symmetric models (NP- 
LR), together with the sensitivity lines for both ATLAS and 
LHCb. 



3 The B s —>J/ifr<f> decay channel 

The B s — > J/*/f0 mode has a clean experimental signa- 
ture. For its reconstruction, the same approach as for the 
Bs,d — > J/i/"7 was used. Details are given in [|6] 

The B s — » J/iff(f> leads to three final state helicity config- 
urations and their linear combinations are CP eigenstates 
with different CP parities [0. For the extraction of the 
CP- violating weak mixing phase cf> s = arg(V* s V c b/V cs V* b ) 
the helicity amplitudes need to be separated. Experimen- 
tally, one can measure three independent angles and the B s 
proper decay time. The initial B s flavour can be tagged in 
part of the events. The ATLAS precision for these mea- 
surements was determined by detector response simula- 
tions and was used as input to angular analyses based on 
a maximum likelihood fit in repeated Monte Carlo experi- 
ments. The likelihood function X. is defined as: 

~ 'J £L t Wit-V) ■ Res(t',t) dt' dt 
with: 

'Wit, Q) = e x co + it, Q) + exofit, Q) + be^ ' (1) 

where e\ - 62=0.5 for untagged events, e^l-r, e^-r if the 
B s is tagged as a particle and ei=r, €2=l-r if the B s tagged 
as an anti-particle, r is the wrong tag fraction, b is the level 



of background and To is the average decay rate of the back- 
ground. The probability density oj* is defined by: 

in which f(t) i +(_) are expressed in terms of four param- 
eters of B s weak decay and oscillations: AT S , T, x s = 
Am,r, and four independent parameters of helicity am- 
plitudes: \A\\it = 0)|, \A±it = 0)|, and the two strong phases 
S\ and 8%. The angular functions Fj are given in f 1 1 01 1 1 11 - 
where more details about the analysis are also given. In the 
maximum likelihood fit the Ar s , the T and the weak phase 
(p s were simultaneously determined along with two helic- 
ity amplitude values and their strong phases. The mixing 
parameter x s was assumed to have been measured in the 
B s — > D s 7r and B s — > D s ai and was fixed. 

While all the eight parameters are independent in the the- 
oretical models, the experimental resolution causes some 
of them become correlated, particularly the two strong 
phases, which prevented their simultaneous determination. 
Thus in the final analysis these two parametres were fixed. 
For an integral luminosity of 30 ftr 1 , AT, can be detrmined 
with a relative error of 12%, while the precision of <p s de- 
pends on the value of x s and on the proper time resolution, 
as shown in Fig. 3, where the discovery lines for ATLAS 
and LHCb are displayed in the (x s - <f> s ) plane, together 
with regions allowed by the Standard Model and by one 
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example of new physics known as the Left-Right symmet- 
ric model r i!2l . 

4 Conclusions 

ATLAS is expected to measure several parameters in the 
Bs meson system, which will be valuable input for the 
flavour dynamics in the Standard Model and its possible 
extensions. 

A clear signal of the decay mode B s — > J/ifrrj is expected 
to be reconstructed with the ATLAS detector, and several 
measurements are possible such as the measurement of the 
branching fraction and the Bs meson lifetime. The exper- 
imental resolution of the CP asymmetry measured shows 
some sensitivity to very large asymmetries as predicted by 
some models beyond the Standard Model. The measure- 
ment of the y angle seems, however to be compromised 
because the Bd — » J/if/r] is not observable. After 3 years at 
luminosity 10 33 cm _2 s _I ATLAS is expected to determine 
AT s with the relative errors of 12%. The ATLAS discovery 
line (Fig. 3) shows that the precision is high enough to be 
sensitive to new physics. 
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